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Abstract. A method is developed to establish the J-resistance function for an isotropic 
linear viscoelastic solid of constant Poisson’s ratio using the single-specimen technique 
with constant-rate test data. The method is based on the fact that, for a test specimen of 
fixed crack size under constant rate, the initiation J-integral may be established from the 
crack size itself, the actual external load and load-point displacement at growth initiation, 
and the relaxation modulus of the viscoelastic solid, without knowledge of the complete 
test record. Since crack size alone, of the required data, would be unknown at each point 
of the load-vs-load-point displacement curve of a single-specimen test, an expression is 
derived to estimate it. With it, the physical J-integral at each point of the test record may 
be established. Because of its basis on single-specimen testing, not only does the method 
not require the use of multiple specimens with differing initial crack sizes, but avoids the 
need for tracking crack growth as well. 


INTRODUCTION 

If a master fracture-resistance function exists for a viscoelastic material, it would -most 
likely- depend on temperature and loading rate. If, in addition, the material has tearing 
ability, that resistance would depend on the amount of crack growth. In either case, any 
reasonably complete characterization of the fracture property function would require a 
large amount of data. Under such conditions, obviating costs becomes important. Such 
is the aim of the work in this document. 

J-resistance, or more simply, J-R curves, characterize the fracture toughness property of 
materials that fail by ductile tearing. In principle, the J-R curve of a material, whether 
elastic-plastic or viscoelastic, may be constructed by testing multiple specimens having 
different initial crack sizes; the curve sought being the locus of collected pairs of initial 
crack size and J-integral at crack-growth initiation. To obviate cost, however, most J-R 
curve testing today utilizes the single specimen technique, which requires knowledge of 
the crack size corresponding to each point of the test record. Traditionally, crack sizes 
have been estimated either by unloading compliance measurements, or by the change in 
elastic potential with crack growth [3]. 

For years, researchers have been developing normalization procedures to increase the 
accuracy of crack-size estimates in the single-specimen technique, especially for tests at 
either elevated temperatures or high mechanical rates. Normalization has its origins in 
the work of Ernst et al [5], who showed that load-separability allows the J-R curve to be 
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estimated using the single-specimen technique with the load-displacement data of a test 
article of known initial crack size. 

In its present form, normalization utilizes the actual load-vs-load point displacement 
record, the initial and final crack sizes, and an assumed analytical relation between 
crack size and displacement to estimate intermediate crack sizes. The method is being 
considered for adoption by several standards. In particular, a draft appendix has been 
prepared and reviewed for standard ASTM E-1820 [2]. 

Although a bounding scheme has been used recently to apply the normalization process 
to polymers and dynamically loaded elastic solids [6], a general approach to treat rate- 
dependent materials is still missing. The work here is concerned with linear isotropic 
viscoelastic solids with a constant Poisson’s ratio. For these materials, an approach is 
developed to construct the J-R function from the single-specimen technique using the 
results of constant-rate tests. In essence, a relation is derived to estimate the crack size 
that corresponds to each load-displacement pair of the test record, and the means to use 
such information to ascertain the corresponding J-integral. 

The method developed here uses the fact that, under constant-rate testing, the J-integral 
at crack-growth initiation for a fracture specimen of fixed crack size may be established 
from knowledge of the crack size itself, the values of load and load-point displacement 
at initiation, and the relaxation modulus of the viscoelastic material, without reference 
to the actual test data at that crack size. As suggested in Figure 1, the method defines a 
relationship that links a generic point on the load-displacement curve of a single test 
specimen with a growing crack to the curve that would correspond to an article of the 
same viscoelastic material but with a crack of fixed size, and to the secant path another 
test specimen of a fictitious linear elastic solid with the same crack size would follow. 
Using the approach of Ernst and co-workers, as well as load and J-integral separability 
of viscoelastic effects for the present class of materials, an equation is derived whose 
solution yields the crack size at each load-displacement pair. 



Figure 1. Physical and linear elastic fictitious paths 


LOAD DEFLECTION RELATIONSHIP IN CONSTANT-RATE TESTING 


The load-vs-load-point displacement response for a fracture test of a linear viscoelastic 
material of constant Poisson’s ratio may be expressed as the product of two separate 
functions. One of these functions accounts for geometry and Poisson’s ratio, and the 
other, for material heredity, in the following form [7]: 

P{a,t) = g(a,L,v) ■ [ E(!—^-)~dT (1) 

o- a r 

In this expression, P and A represent, respectively, the externally applied load and load- 
point displacement, and the symbol L is used to denote all parameters characterizing the 
geometry of the test article, other than crack size, a. Also, E, is the tensile relaxation 
modulus of the material, and v, its -constant- Poisson’s ratio. 

If the fracture test is performed at a constant rate the load-vs-displacement relationship 
resembles that for a linear elastic material. This form is arrived at by direct integration 
of (1), after dividing and multiplying the result by t. Mathematically: 

P(a,t) = g(a,L,v) ■ E(t) ■ A(/) (2) 

Which indeed has a linear-elastic character, except that E represents the time-averaged 
value of the relaxation modulus. That is: 


i t 

E(t) = -- ^E(s)ds 

t r\ 


(3) 


J-INTEGRAL FOR A STATIONARY CRACK IN CONSTANT RATE TESTING 

Reference [7] shows the J-integral at fixed crack size for a viscoelastic solid of constant 
Poisson’s ratio is proportional to the area under the load-vs-load displacement record of 
the fracture test. As it turns out, the factor of proportionality depends on crack size and 
other relevant geometric features of the test article, and perhaps also on the —constant— 
Poisson’s ratio of the material at hand, but not on its relaxation modulus. Thus: 


J = 





(4) 


This expression is similar to that in current use to estimate the plastic J-integral [4]. In it, 
the factor T](a, v,L), is merely a crack-configuration function. It is related to the geometry 
function, g(a,v,L), for a specimen of uniform thickness. Bn, in the following manner [7]: 


( 5 ) 


rj(a,v,L) = - 


a[ln(g>] 

da 


For tests carried out under constant rate, the J-integral may be expressed as the product 
of two functions. One function depends solely on the heredity of the viscoelastic solid, 
and the other is the J-integral that would be obtained for a linear elastic material at the 
current levels of load and displacement. The latter integral is simply called the secant 
J-integral, 

• 

In fact, taking a constant-rate loading history A = A- / and the load-vs-displacement 
relation (2), into (4), the following is obtained, after changing integration variable and 
limits, and pulling out of the integral quantities independent of the integration variable. 


J = 


7}(a,v,L ) 


B 
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■ g(a,v,L) ■ A- A- Jr • E(r) dt 
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Multiplying and dividing the previous expression by t 2 ■ E(t ) and re-grouping: 


j = U ^ R . S^L £t £(,) - A • A ■ - f T • E{ r) dr 
B n t ■ E(t) t J 0 

And, upon multiplying and dividing by 2, and again recalling expression (2): 


J = 


1 f — 

-Jt- E{r) dr 
2- 1 0 




B k , 2 


(6) 


Or, more succinctly: 

J(a,t) = a(t)- J s (a,t) ( 7 ) 

Where the function «(/), which accounts for material viscoelasticity, and the secant J- 
integral, J s , corresponding to a —fictitious— linear elastic solid, have been defined as: 

- \r E(r)dr 

a(t) = 2-— — = ( g ) 

t-E(t) 


and: 
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(9) 


J s (a,t) = 


2 B k 


r}(a,v,L)- P(a,t)- A(t) 


The form in (7) is reminiscent -and equivalent- to the relationship between the physical 
J-integral, and the viscoelastic J-integral, J v , introduced by Schapery to characterize the 
initiation fracture toughness of viscoelastic materials [8], Because the secant J-integral 
is expressed directly in terms of physical quantities, it is easier to work with than J Y . 


J-INTEGRAL FOR A GROWING CRACK IN CONSTANT RATE TESTING 

By definition, the values of P and A are known at every point of the load-displacement 
record, irrespective of whether or not the crack size, a, is known. On the other hand, 
the crack-configuration function, 7} , is known -and constant-only up to the onset of 
crack-growth initiation. From there on, it is not known a-priori what crack size, a, is 
associated with any given test pair (A,P); a fact that prevents direct use of expression 
(9) to estimate J s , and hence, J. This difficulty may be circumvented by taking the total 
differential of (9), considered as a function of a, and t, and by integrating the resulting 
expression along paths of constant a, and constant t, as shown subsequently. 


From (9), the total differential of J s becomes: 
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dJ (a, t) - 
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dr} dP 

— -PA + 7} —A 
da da 
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dt 


The indicated partial derivatives, evaluated using (2), are: 
d S ~B A = I4g „.r. 


dP_ 

da 


E- As • g ■ E ■ A = -7} (a, L, V) ■ P{aj)[ 

da g da 


and: 


dP 

dt 


d_ 

dt 


g E- At 


= g ■ A~jE(s)ds = g-A-E(t) 


Combining these expressions and recalling the load-deflection relationship (2), the 
definition of J s , and that constant-rate testing is involved, the following results: 


dJ s {a,t) - y(a,v,L) ■ J s (a,t)\da + 
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t-E(t) + j E(s)ds 


dt 


(10) 


In which the following definition has been introduced for convenience: 



( 11 ) 


y(a,v,L) 


1 drj_ 
rj{a,v,L) da 


7](a,v, L) 


Expression (10) may be used to obtain J s at any point on the load-displacement record 
using any convenient load path. In particular, since under constant-rate testing 1 and A 
are proportional, any generic point (A ,P), or ( t,P ) of the test record may be reached, as 
suggested in Figure 2, by integrating along a path of constant crack size, followed by 
one at fixed time -or fixed displacement. Indeed, reaching point (t,P) along path OBC 
in the figure: 


it 

J s (a,t) = jy(a,v,L)J s (a,t)\da 


+ 



j t- E(r)dr 
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N Ji:(.s)ds \dz 
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Figure 2. Fictitious integration paths to adjacent points of the test record 


Whose second term may be cast in a more useful form. Integrating by parts the first 
integral inside the brackets, collecting like terms, and, using the load-displacement 
relation (2), multiplying and dividing the result by t, the whole second term becomes: 


l = J s (a„t ) = %^/>(a„f)-A(r) 


2 • B, 


And hence: 


a 

J s (a,t) = \y(a,v,L)J s {a,t)\da + J 5 (a,,t) 

a , 


Or, equivalently, because time and displacement are proportional: 


J s (a,A) = jy(a,v,L)J s (a,A)\^da + J s (a,, A) 

a, 

So that, if test point (4+;, P,+i) is reached along path a = a„ followed by path A = A i+t : 


u <+ 1 

J s (a M A M )= \y(a,v,L) J s (a, A M )da + J s (a,, A, +1 ) 


(12) 


Since along a path of constant crack size the slope of the load-deflection line is P/A: 
JM .A,.,) = ^A-P(a, A W ).A,„ . A-(%) 
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Which, upon multiplying and dividing by (4)" becomes: 


/ A \ 


J s (a,A M ) = J.MA) 

And, similarly: 
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(13) 


As it should be, since the secant J-integral corresponds to a linear elastic solid, however 
fictitious. 

With the previous results, (12) may be expressed more conveniently as: 


M i+ 1 

7>,. +1 ,A, +1 ) = J y(a) • 7 ,( 0 , A (+] ) • da + 






J s (°n A,) 


(14) 


Note that, although a, is known, a i+t is not. However, if 4 is selected close to 4+ 1 , a, and 
a (+/ will be close together. For this reason, the integral in the previous expression may be 
evaluated approximately. Using the trapezoidal rule for this purpose, invoking (13) to 
convert J f (a, ,4 +/) to /*(«, ,4), and (9) to express the resulting secant J-integrals in terms 
of their corresponding rj functions, the following results after rearrangement: 


V(a M ) = 
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Since the form of the crack configuration function, 77 (a), may be established by analysis 
of the test specimen -either exact or numerical- the function }(a) may be obtained from 
it through (11). With this, expression (15) may be solved iteratively for the crack size, 
a i+ 1 , corresponding to each test point (P i+1 , 4 +/)- This requires the triplet (P<» A a , a () ) to 
be available at the outset. Incidentally, because a, is close to a (+/ , the former constitutes 
a good initial estimate of the latter in the iteration process. Knowledge of the crack size 
at each point of the test record permits the secant J-integral to be estimated per (9), and 
the corresponding physical J-integral via (7) through (9). 

The set of triplets, (a„ t„ /,), associated with each constant-rate test will define a portion 
of the J-resistance surface, J(a,t), similar to that shown with the solid line in Figure 3. 
The collection of triplets from a sufficient number of tests carried out at different, but 
constant, rates will define the J-resistance surface. Clearly, if the viscoelastic material in 
question is thermo-rheologically simple, the time-temperature superposition may be used 
to stretch the time scale of the tests in the usual manner [ 1 ]. 



Figure 3. J-Resistance Surface of a Viscoelastic Solid 


CONCLUSION 

The present work is concerned with isotropic linear viscoelastic solids having a constant 
Poisson’s ratio. For such materials an approach is developed to construct the fracture 
resistance surface using the single-specimen technique with constant-rate tests. The basis 
of the proposed method is that at constant test rate and Poisson’s ratio, the J-integral at 
crack-growth initiation for a specimen of fixed crack size is proportional to the J-integral 


of a fictitious linear elastic solid at the same values of crack size, load, and displacement. 
In other words, it is not necessary to have knowledge of the complete load-vs-deflection 
curve at each crack size to obtain the corresponding initiation J-integral. The fictitious J- 
integral may be obtained once the crack size is known at each point on the curve. Using 
that the loading paths corresponding to the fictitious J-integrals are secant lines drawn 
from the origin, an equation is developed linking the crack size at a generic point to that 
of a previous point close-by. Solution of this expression yields the crack size, and hence 
the fictitious J-integral, and from it, the physical J-integral sought. The approach may be 
repeated for as many points of the test as are necessary to resolve the J-resistance curve. 

The method does not require use of multiple specimens with differing initial crack sizes, 
and avoids the need for tracking crack growth as well. The constancy of Poisson’s ratio 
may be relaxed using specimens whose response is independent of it —such as single- 
edge-notch beams- or depends on it only weakly -as with compact tension plates. 
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